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Abstract. Quantum effects in weakly disordered systems are governed by 
the properties of the elementary interaction between propagating particles and 
impurities. Long range mesoscopic effects due to multiple scattering are derived by 
iterating the single scattering vertex, which has to be appropriately diagonalized. 
In the present contribution, we present a systematic and detailed diagonalisation 
of the diffuson and cooperon vertices responsible for weak localization effects. We 
obtain general expressions for eigenvalues and projectors onto eigenmodes, for 
any spin and arbitrary elementary interaction with impurities. This description 
provides a common frame for a unified theory of mesoscopic spin physics for 
electrons, photons, and other quantum particles. We treat in detail the case 
of spin-flip scattering of electrons by freely orientable magnetic impurities and 
briefly review the case of photon scattering from degenerate dipole transitions in 
cold atomic gases. 



PACS numbers: 73.20.Fz, 03.65.Fd, 72.10.Fk 

1. Introduction 

The physics of multiple scattering is governed by the iteration of elementary scattering 
events. The description of mesoscopic effects, due to phase coherent multiple scattering 
of waves, therefore requires the elementary interaction to be in a form which is 
suitable for the iteration. The scattered particles in question may be electrons, 
photons, neutrons, or cold atoms; the scatterers may be point-like impurities, spin-flip 
impurities interacting with the electron spin via an exchange interaction, spin-orbit 
impurities, atoms interacting with the photons via the dipolar interaction, or classical 
dielectric light scatterers £Q , to name a few. 

On a classical level, multiple scattering is described by a Boltzmann-type 
transport equation, which in a microscopic description is generated by considering 
pairs of complex conjugate amplitudes co-propagating along the same scattering path. 
In a diagrammatic representation, these amplitudes are depicted by the so-called 
ladder diagrams. The sum of these ladder diagrams constitutes the "diffuson" which, 
in the long-distance limit, describes a diffusion process. Weak localization corrections 
to classical diffusive transport are described by the "cooperon" , the sum of so-called 
maximally crossed diagrams made of amplitudes that are counter-propagating along 
the same scattering path. If the wave scatters off structureless point scatterers, the 
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maximally crossed diagrams can be disentangled by returning one amplitude and 
thus transform into a sum of ladder diagrams. Then, reciprocity |2] assures that the 
quantum correction due to this interference is maximal. 

In the presence of scatterers with internal degrees of freedom, multiple scattering 
contains richer physics, since these internal degrees of freedom couple to the degrees 
of freedom of the propagating wave. The two subsytems, propagating wave and 
impurities, become entangled, and discarding all which-path information by tracing 
out the unobservable impurity degrees of freedom leads to an effective dephasing of 
coherent effects for the observed wave. In mesoscopic electronic samples, for example, 
the spin of a propagating electron couples to the spin of magnetic impurities, and in 
light-scattering atomic clouds the polarization of propagating photons interacts with 
the internal atomic angular momentum. In these cases, the elementary scattering 
vertex in the diffuson series is a tensor with 4 spin indices, connecting two incident 
spin states to two scattered spin states. A succesful derivation of multiple scattering 
properties then requires that this elementary vertex be iterated. This problem has 
been studied and solved in several specific cases. For instance, in the case of spin- 
orbit coupling and scattering by magnetic impurities, Hikami, Larkin and Nagaoka 
showed that the cooperon can be diagonalized in the singlet and triplet subspaces. 
Similar methods were employed for calculating conductance fluctuations 0]. In the 
context of light scattering by thermal atomic gases, the iteration structure of the 
photonic diffuson was completely determined by Barrat, Omont and others |SJ|S1|7]. 
For the study of phase-coherent effects in cold atomic gases, the cooperon for the atom- 
photon problem was calculated exactly by two of us [Hj- Its diagonal properties were 
used to describe coherent backscattering by cold atomic gases and weak localization 
phase coherence times {5]. In all these cases, the iteration of the elementary vertex 
properties was done by hand, finding the appropriate diagonal tensors and associated 
eigenvalues rather heuristically. 

In the present contribution, we provide a thorough understanding of the diffuson 
and cooperon vertex diagonalization, of the different projectors and eigenvalues 
involved, of their spin and coupling dependence as well as of their precise relationships. 
In the next section we start by recalling the example of the spin-flip scattering of spin 
^ particles. In section |3J we present a general diagonalization scheme for arbitrary 
ladder and crossed vertices. We derive in detail the isotropic projectors onto invariant 
subspaces for scalar vertices. Once the algebraic structure has thus been laid, we 
calculate in section 0] the correponding eigenvalues from the microscopic scattering 
potential. Finally, we conclude this paper by indicating some possible extensions of 
the work. |Appcndix A| contains a brief review of photon scattering properties in the 
light of the present work. 

2. A heuristic diagonalization of the electronic spin-flip vertex 

2.1. Diffuson and cooperon 

Let us first consider multiple scattering of a quantum particle by randomly 
distributed impurities without internal structure [Q. The disorder-averaged 
probability of a wavepacket emanating at point r and time t = with 
density matrix po( r ) to be detected at point r' and time t > is given by 
P(r,r',t) = {(r'\U(t)po(r)W(t)\r')} where (. . .) av signifies a trace over the impurity 
configurations. For a quasimonochromatic wave packet of central energy e and long 
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Figure 1. Three equivalent representations of the cooperon interference at five 
scatterers: (a) Real-space representation of counter-propagating amplitudes; (b) 
momentum-space maximally crossed diagram with retarded propagator (upper 
full line) and advanced propagator (lower dashed line) connected by impurity 
scattering events (dotted lines); (c) the same diagram with returned advanced 
propagator line exhibiting the ladder structure. 



evolution time, the Fourier transform of the detection probability is P(r,r',u>) = 
[2Trp(e)]~ 1 (G R (r, r', e)G A (r, r', e — a;)) in terms of the retarded and advanced 
Green functions G R ' A (E) and the average density of states p{e). 

This average probability satisfies an integral equation of the Bcthc-Salpeter 
type generating a multiple scattering sequence. In weakly disordered samples, the 
interference of amplitudes propagating along different scattering paths will be washed 
out by the disorder average. Therefore, the dominant contribution will come from co- 
propagating amplitudes along identical scattering paths, thus discarding interference 
effects and recovering a classical propagation picture. This propagation is described 
by the so-called diffuson, the propagation kernel of the multiple scattering sequences 
defined in operator form by D = L + LGD where L describes the elementary scattering 
by a single impurity, and the four-point operator G = (G A ) av (G R ) av is the intensity 
propagator (of Boltzmann-Drude type with factorized averages) between scattering 
events. In diagrammatic representations, this series has a ladder structure and 
can formally be summed as a geometric series, D = L/(l — GL). Going to the 
diffusion approximation (Kubo-limit of large distances and long times) permits to 
derive the effective diffusion constant of this classical diffusion process as function 
of the microscopic parameters, in a spirit similiar to the kinetic equation in the 
Boltzmann-Lorentz model of classical particles colliding with fixed impurities. In 
the case of electrons, the Einstein relation between the diffusion constant and the 
conductivity then allows to recover the classical Drude conductivity. 

Quantum corrections generated by the interference of amplitudes propagating 
along different scattering paths can be incorporated by considering more general 
scattering processes. A particular interference that survives the disorder average is 
the weak localization correction of the classical diffusion constant due to amplitudes 
counter-propagating along identical scattering paths ^Oj) depicted in figure^a). In a 
diagrammatic represention, this interference is given by maximally crossed diagrams 
[figure njb)] that can be unfolded to a ladder structure [figure [TJc)] and summed to 
C = X/(l — GX) where X is the single scattering vertex for the returned advanced 
amplitude line. For simple impurities without internal structure, one has X = L. 
Weak localization then enhances the classical return probability of a particle by a 
factor of 2 and reduces the diffusion constant. Experimentally, this effect can be 
measured for instance in the electronic negative magnetoresistance where an external 
magnetic field suppresses the weak localization corrections and thus leads to a larger 
conductance. In optics, the interferential enhancement of backscattered intensity is 
called coherent backscattering and has been observed in a large variety of samples. 
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The case of impurities with internal structure acting on the spin degrees of freedom of 
the propagating particle is more involved as should become apparent in the following 
example of electronic spin-flip scattering. 



2. 2. Definition of the spin-flip vertex 

We consider now a particle of spin S propagating in a disordered sample with spin J 
magnetic impurities. The particle spin states are written as eigenstates \sa) of S 2 and 
S z . For electrons, s = ^ and a = + \. The spin operator then is S = er/2 (in natural 
units h = 1) where the components of er = (a x , <r y , a z ) are the usual Pauli matrices. 
The interaction with a given magnetic impurity is described by the hermitian operator 
V m — g J ■ S with coupling strength g. The Born scattering amplitude for the spin-flip 
process \sa) t— > \sj) is 

(sj\V m \sa) = gJ ■ (s-f\S\sa) = g J ■ S ya . (1) 

The effective scattering intensity of this process is described by the four-point vertex 
(the overline indicates complex conjugation) 

C aP , lS = TP= U S1 \V m \sa)( s 5\V m \sf3)) = |g| ' J( / + 1} S ia ■ S ps (2) 

P-Li-5 X /av 3 

where (. . .) denotes a trace over the impurity configurations, here an isotropic 
average {JiJj) av — SijJ(J + l)/3 over all possible orientations of the freely orientable 
magnetic impurity. This average is the fundamental reason for the non deterministic 
dcphasing of the multiple scattering process. We choose to normalise the spin-flip 
scattering strength g such that the intensity vertex is written as 

5 7 q ■ S/3S . , 

Ca ^ 5 - 8(8+1) ■ (3) 

This normalisation choice endows the vertex with the convenient trace-preserving 
property C aa ,js — here as in the following, summation over repeated spin indices 
is understood. 

Weak localization corrections to transport are embodied in the so-called cooperon 
and are generated by maximally crossed diagrams. These diagrams can be "unfolded" 
to a ladder structure such that the intensity vertex (J3J| is replaced by the crossed 
vertex 

X a f3,~/8 = £ a s,~ff3 = Y = (sjs8\X\sasf3) = 7 , a S f . (4) 
(3-ll-S s{s + l) 

For reasons that will become clear later, the crossed vertex X is denoted by a roman 

letter, whereas the vertex C is held in curly script. 



2.3. Elementary diagonalization 

The natural coupling scheme for the vertices 10) and (@J is the "vertical" combination 
(aj) <-* {(36) between the elementary scattering amplitudes. However, in multiple 
scattering diagrams, the above intensity vertices have to be chained "horizontally" 
in the direction (af3) <-> (7J) according to the following product rule of four-point 
vertices: 

(GC)ap.jS = Cap^v^nv^S (5) 

In the product definition, the order of operators is inverted since operators are 
conventionally applied to their arguments from the left, but their vertex symbols 
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are usually added to a diagram on the right. The rank-four tensor G a p n s = 

(G„ 7 ) av (^Gps"j describes the average propagation between scattering events. By 

virtue of rotational invariance, the average propagators (G Q7 ) av = (G) av (5 Q7 are 
proportional to the identity in spin space such that G a /3 tl s — (C R ) av (G A ) av 8 a ^8fjs 
is proportional to the "horizontal" identity. (Note that for photons, however, 
transversality implies that G a p n s is not the identity which leads to the more 
complicated scenario described in appendix A, featuring nonetheless the general 
properties discussed in the present section.) But in order to calculate the summed 
diffuson V = C/(I — G£) and the cooperon C = X/ (1 — GX), the vertices C a p^s and 
have to be diagonalized with respect to the horizontal direction (a(3) <-* (■yS). 
For electrons with s — i, this amounts to diagonalising 4x4 matrices [I.. It turns 
out that the diffuson and cooperon vertices can be cast in the form 

£ a f3 n 8 = A T a °0 n s + A l ^af^S ' ( 6 ) 
X a f3,-f8 = XO T afj n s + Xl T a0,y6- ( 7 ) 

Here, the diffuson vertex tensors 

1 



TalnS = 7,5^6, (8) 



Tafias — 2 <T ' 5 " ' a '< S ~ SyaSfJS — ~ S/3 a S^S (9) 

are orthogonal projectors with respect to the horizontal product rule 

^Slriv'i = S KK' T «f,\s- (10) 

Likewise, the cooperon vertex tensors 

T alyS = \ ( S l*S Sl 3 -Ssa^p), (H) 



' a/3, 7<5 2 

T aPni = 2 ( 5 ~< a6s P + S Sa^fi) (12) 



are orthogonal projectors such that T^^ U T^ S = S KK >T^ yS . Both sets of 
projectors sum up to the identity 8 ia 8fjs for the horizontal product rule (JSJ. Obviously, 
the diffuson projectors (JSJ and (J5J) are different from the cooperon projectors (fill) and 
H12fl . This is in sharp contrast to the case of photon scattering (s = 1) by atoms with 
degenerate dipole transitions, where the same set of orthogonal projectors can be used 
for both vertex types |S] (see | Appendix Aj for details). 

In the diagonal decomposition JSJ , the eigenvalues of the diffuson vertex are found 
to be Ao = 1 (non-degenerate) and Ai = — | (3- fold degenerate). The eigenvalues 
of the normalised crossed vertex are xa = — 1 (non-degenerate) and Xl = \ (3- 
fold degenerate). It has been noticed that the cooperon spin vertex eigenvalues xk 
correspond to the singlet channel K — and the triplet channel K = 1, respectively, 
which accounts for the degeneracies (3J- Remarkably, the eigenvalues Xk and xk are 
equal in magnitude but opposite in sign, which is not properly explained on this level 
of heuristic diagonalization. 

Prompted by these observations, we wish to answer the following questions: 

(i) Given diffuson and cooperon vertices for arbitrary spin s, which are the orthogonal 
projectors that assure a least redundant diagonalization? 

(ii) How do the diffuson and cooperon eigenvalues depend on the microscopic spin 
scattering mechanism? 
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2-4- First answers 

2.4-1- General idea In essence, the intensity vertices map two incident spins s 
onto two final spins. Furthermore, they are scalar objects since they are obtained 
by an isotropic average over microscopic degrees of freedom. The invariance under 
rotations is then responsible for the eigenvalue degeneracies. The key idea is to 
decompose the argument and image spaces into irreducible subspaces with respect 
to the rotation group. The relevant subspaces are labelled by the effective recoupled 
spin K = 0, ... ,2s. A generic scalar vertex A can only connect irreducible subspaces 
with equal K, and its eigenvalues are degenerate in each subspace. In the appropriate 
recoupled basis, a scalar vertex takes the diagonal form 

( aolo ... ^ 
aili ... 



A 



V ... a 2s l 2s j 



2s 

\ " „__m(K 



a K 



T (K \ (13) 



K=0 



The projectors T^ K ' are simply projectors onto the irreducible subspaces. Therefore, 
rotational symmetry alone dictates that there are at most (2s+l) different eigenvalues, 
each (2K + l)-fold degenerate. This decomposition is optimal if the scalarity is the 
only information available and holds for arbitrary spin. 

2-4-2. Recoupling schemes Clearly, the natural coupling between spin indices in the 
ladder and crossed vertices and Q is the "vertical" coupling scheme (0:7) «-> ((36) 
that is inherited from the scattering amplitude Q). Unfortunately, this coupling is 
not suited for an iteration with the product JSJ. For the diagonalisation, we therefore 
have to recouple the spin indices: into the "horizontal" coupling scheme (a(3) <-* (7J) 
for the ladder vertex, and the "diagonal" coupling scheme (a6) *-* ((36) for the crossed 
vertex. Consequently, all vertex eigenvalues we derive will feature 6j-symbols that 
describe the recoupling of 4 spins in angular momentum theory. 

By exchanging the indices 6 <-> (3 in the crossed vertex, we are actually able to 
recover formally the ladder structure, but there is a price to be paid. Taking seriously 
the disposition of kets and bras, we see that the two vertices have different rotational 
structure: the crossed vertex defines a linear mapping between product states 

X : \sa)\s(3) -> \sj)\s6) , 

whereas the diffuson vertex is a mapping not between states, but between operators: 

C : \sa)(s(3\ h-> \sj)(s6\ . 

This difference, due to the exchange (s(3\ <-> \s6) and subsequent relabeling (3 *-* 6 
in the unfolding procedure of the crossed to the ladder series, leads eventually to 
two distinct sets of projectors. We will therefore diagonalise the diffuson vertex as a 
superoperator, but treat the cooperon vertex as an ordinary operator. 

2.4-3. Projectors In sectional we will explicitly construct the diffuson projectors onto 
irreducible subspaces. For spin |, this reduces indeed to the electronic projectors © 
and ©■ The cooperon vertex projectors l|ll|) and (|12fl will be shown to be given by 

T (o) _ _ i T (o) l r (i) ru x 



r (i) = 3 T (o) l r (i) 

~otj3,"f8 2 aS,-yj3 ' 2 Q^,7/5 



T (i) _ 3 T (o) , l r (i) 



Mesoscopic scattering of spin s particles 7 
This is in fact the spin i version of the more general relation 

K 1 

between the cooperon and diffuson projectors, valid for arbitrary spin, to be derived 
in section ?6. 61 Here, our notation 

R S (K,K') = (2K+1)( I S s * | (17) 

is a 6j-symbol from standard angular momentum theory [111 112) . Thanks to the 
6j-symbol orthogonality (35c)] $ 

^R.(K,K')R,(K',K")=8 k ,k», (18) 

K' 

the inverse relation to l|16fl is equally simple: 

K' 

2.4-4- Eigenvalues We will show in section 0] that the eigenvalues of the normalized 
spin-flip vertices (0) and Q are given by 

The eigenvalues for an arbitrary microscopic spin scattering potential will be derived 
below in full generality. In all cases, the eigenvalues of a scalar scattering vertex are 
linked by the recoupling relations 

\ K =J2 R s(K',K) X k>, xk = Y j R ^ K "' K ) X K"- ( 21 ) 

K' K" 

These relations between eigenvalues take the form of a contravariant transformation 
of coordinates associated with the respective covariant transformation and (|19H 
of the projectors. 

These results should motivate our readers to consider with interest the following, 
more involved derivations. In the following scction|31 we lay the algebraic foundations 
of the decomposition by deriving the orthogonal projectors, before turning to the 
eigenvalues in section 0] 

3. Diagonalization of intensity vertices 

Let us now consider a general spin interaction defined by its matrix elements (s7| V|sa) 
for arbitrary spin s. The corresponding diffuson and cooperon vertices are given by 

C a p t75 = (( S1 \V\sa)(sp\Vi\s5)) av , (22) 

« S7 |^| S a)< S( 5|^t| s /3)) av . (23) 

X Instead of compiling a large appendix, we will refer to standard definitions and sum rules by citing 
the exact location in appendix C of Messiah's book, e.g., his equation (35a) by writing 1121 (35a)]. 
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3.1. The diffuson vertex as a spin superoperator 

Spin states \sa) are vectors in the Hubert space H s = <C ds with dimension d s = 2s + l. 
The diffuson vertex £ is a linear mapping |sa)(s7| i— * \sS}(s/3\ between spin operators. 
Its argument space therefore is the space of linear operators acting on H s , the so- 
called Liouville space L(H S ) with dimension d 2 s Any linear operator A E L(H S ) 
is simply a d s x d s matrix. The trace-preserving vertex £ : L(H S ) — > L(H S ) then 
is a superoperator (thus the notation with a curly script), mapping a matrix A onto 
another matrix A' = LA. Its action in the basis of spin projectors {|sa)(s/3|} reads 
A' s — C a p n sA a j3 in terms of the matrix elements 

£^ )7 5 = tr{(|^)( S7 |)£(|sa)(^|)} . (24) 

Here tr{-} = ^2 a (sa\ -\sa) is the trace over H s . In superoperator notation, the diffuson 
vertex reads 

C = E (\<n)(* s \)£*W tr{(| S /3)( S a|) • } (25) 

af3jS 

With this notation, the ressemblance with the Liouvillian £ = — i [H, ■ ] , the generator 
of time evolution, becomes apparent. We define the trace of £ as a linear operator in 
the Liouville space as 

Tr L £ = ^ £af3, a/3 • (26) 

a,/3 

In the extensive literature on Liouville space formalism (141 1151 llfil 117) . one 

often views the spin operators as vectors in L(H S ) and defines corresponding kets 
by |s/3)(sa| = \aj3)). Using this notation, the diffuson vertex matrix elements are 
given by C a p,-yS = {{l5\C\aj3)) , and the superoperator takes the very simple form 
£ = J2 a 0jS ^q/3,7<5|7^)) (( a P\- I n this notation, the parallel with the cooperon vertex 
operator (see Q44JI below) is especially clear. However, we deliberately choose to use 
the superoperator formulation in the following because it allows us in section 13.3.21 
to derive the ladder vertex projectors in terms of spin operators (which is needed to 
get expression (0). Moreover, to answer completely question 1 raised in section I2"31 
we have to explain the difference between the diffuson and cooperon eigenstructures 
rather than their similarity. This difference reflects the different behaviour under 
rotations of spin states \sa) and their conjugates (sa\ that are explicitly featured in 
the superoperator notation 1)25)1 . 

3.2. Decomposition into irreducible superoperators 

An incident state \sa) in scattering amplitudes like is a spinor, i.e., a vector 
in H s that transforms under the irreducible representation of the rotation group 
SU(2): \sa) U\sa) = \sjj,)(s/j\U\sa) = ^ t/ Ma |s/x) with an appropriate unitary 
and unimodular matrix UW — td B , detU = +1. A final state (sj\, however, 
transforms contragradiently [T^], i.e. , under the complex conjugate representation 
D( s \ (s-y\ i— > (s^\U' = U vl (sv\. Under a rotation, the complete spin vertex is 
transformed as Cap^s *—> U a s{U*) lp C lJLUt p Cr U lJLa {U^)i3 v . Clearly, this vertex is not a 
rank four tensor (that would transform under the direct product (2?( s ))® 4 ), but rather 
a two-by-two mixed tensor that transforms under (D^ ® D^)® 2 . If the vertex is a 
scalar, it is invariant under this transformation. 
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In expressions 11240 and l(25|) . the operator arguments of C are decomposed 
over the decoupled product basis {|sa)(s/3|} of the Liouville space L{H S ). But 
according to our diagonalization strategy, we want to use a basis adapted to irreducible 
representations of the rotation group. We first perform the Clebsch-Gordan (CG in 
short) decomposition of the argument and image representations ® DO). These 
are then recoupled in turn to give the complete CG-decomposition of {D^ s > eg) DO))® 2 
which yields the irreducible components of the superoperator. The route thus taken 
may be traced in the following map: 



(£)W ® DO)) (D^ <g> DO)) 

\ / \ / 

J)^ ® D^ K ') (27) 

\ / 

D (L) 

Following standard procedures from angular momentum theory |18II13| . one can define 
a set of irreducible operators adapted to our purpose, 

T (K) = T ( K \ s ,s) = J2(-y~ m ( ssm '- m \ K l) \sm')(sm\ (28) 



in in 1 



with matrix elements 

(sm'\T^\sm) = (-) s - m {ssm'-m\Kq) (29) 

where (ssm! —m\Kq) are the usual CG-coefficients. These types of tensors, called 
"statistical tensors" or "state multipoles" are irreducible components of the density 
matrix and have been developed by U. Fano and G. Racah in the 1950s |18| . 
Their construction is very similar to the coupling scheme of angular momentum 
eigenstates: one simply chooses a linear combination of spin projectors with suitable 
CG-coefficients. The definition (128(1 features a characteristic minus sign in front of the 
spin quantum number m that is reminiscent of the contragradient transformation of 
(sm\. Hermitian conjugation is defined by Tq" = (— ) 9 Tig . 

The orthogonality of CG-coefficients assures that the operators 1(28(1 are 
orthonormalized with respect to the matrix scalar product (A\B) = tr{A^B}, 

(T^\ if >) = tr{T^Tp} = 5 KK ,5qq, . (30) 

The set of irreducible tensor operators Tq provides a natural basis that incorporates 
best the rotational symmetries. Any linear operator O can be developed in this basis 
according to O = J^Kq OKqTq K \ with components 

Kq = (tW| O) = ^{-) s - m { S sm'-m\Kq){ S m'\0\sm). (31) 

mm' 

Inserting \sm')(sm\ — Y^Kq(~ ) s ~ m (ssm'— m\Kq)Tq K ^ in the vertex definition 125(1 . 
the superoperator becomes 



c= V TpL q % ql (Tf\ •) (32) 



KqK'q' 



(K') 



CT ( q K) ) or 



where its left-right irreducible components are C q K ^, q , = 
4% = E H a - p (ssa-0\Kq) C a0n& (-)*-* {ss8 - 7 \K'-q'). (33) 

afi-yS 
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Now we recouple the irreducible argument and image representations to get the 
complete CG-decomposition (last line in 1)27(1 1. We define a basis of irreducible 
superoperators Tm (K, K') of rank L with components m = —L, . . . , L: 

Ti L \K,K<) = J2(-) K - q (K'K q '- q \Lm)Tp (t^\ ■ ) (34) 

qq> 

The recoupled objects can be precisely located on the decomposition map l|27|) : 

\sa) (afi\ \s6) (s 7 | 

tJ*> TP (35) 



\ 



tL l \k,k') 



The vertex in irreducible superoperator notation is then 

£ = EE K '^ L) ( K > K ') ( 36 ) 

KK' Lm 

with components C Lm (K,K>) = E qq >(-) K ~ q ( K ' K <l'-<l\ Lm )£ q % q >- 
3.3. Scalar diffuson vertex 

The above basis set construction and decomposition into irreducible superoperators 
apply to arbitrary superoperators. This basis change is especially profitable when the 
diffuson vertex under consideration is a scalar with respect to rotations. In this case, 
its only non- vanishing irreducible component is Cqo(K, K') for L = 0, m = 0. The 
usual selection rules of CG-coefficients then require in (|36[1 that K — K' and q = q': 
scalar superoperators indeed connect irreducible subspaces L(H S )^ K ^ with equal rank 
K. Each of these subspaces has dimension (2K + 1) and the total dimension is of 
course preserved, J2k=o(^K + 1) = ^s = + I) 2 - As found heuristically in section 
12.31 L then takes the form 

2s 2s 

C = A>o (K, K)T m (K, K)=J2 A^T W (37) 

K=0 K=0 

with X K = C 00 (K,K)/y/2K + 1 and = ^/2K + 1 T (a) (K,K). The calculation 

of eigenvalues will be treated in detail in section 0] Wc now complete the algebraic 
characterization of the projectors. 

3.3.1. Properties of the scalar projectors The superoperators T^ K ^ which diagonalize 
a scalar vertex are projectors onto the Liouville subspaces L{H S )^ K ^ of irreducible 
operators of rank K: 

T W = V2K + 1 T (0) (K) = T [ q K) ( T W | ' ) - E T v K) tr { T P ■ } ( 38 ) 

9 1 

These operators are scalar objects themselves since 

E/f'rf" generalizes the 
scalar product between vector operators (K = 1) to arbitrary rank K (87)]. 
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Thanks to the orthogonality l|3U|> of the basis tensors, the T^'s are indeed orthogonal 
projectors, 

T (K) T (K') = £ T W (TWlT^) (t^I • ) = 6 KK ,TW . (39) 



qq' 



°KK' & qq > 

Their matrix elements in the decoupled basis {|sa)(s7|} of spin projectors are found 
by inserting l|38|) into the superoperator definition l|24fl. 

^-tr{(|5)( 7 |)TW(|a)(/3|)} =^( 7 |Tf)|J)(/3|Tf )t )a) . (40) 

Using the matrix elements 1)291) and the completeness relation of CG-coefficients 12 
(14a)], it is straightforward to show that 

Tr L TW = ^T a ^=2K + l, (41) 

a/3 

as expected for the identity in the subspace L(H S )^ of dimension 2K + 1. 
Furthermore, the projectors T^ K ' sum up to the identity with respect to the horizontal 
product rule ©: J^K^lf^s = Ia/3,-?S = 5 ia 8f jS . 



3.3.2. Expression in terms of spin operators The projector onto scalar operators is 
the "trace-taker" 

r<°) = 1 1(1, • ) = j- 1 tr{ • } , (42) 

with 1 the identity in H s , which is all but a surprise considering that the scalar 
part of a matrix is its trace. In the decoupled basis, we have T^i ' s = -^-8 a p5 1 s- 
justifying thereby the tensor (JSJ found by elementary diagonalization in the electron 
case. Now it is evident that a unit superoperator eigenvalue Ao = 1 is equivalent with 
trace-preservation (and hence particle/energy conservation). 

But already for the projector onto vector operators, using 1)40(1 involves 

a sum over products of CG-coefficients, and it is advisable to look for a more 
transparent formulation. Equation 1(38)1 tells us that we need to find a contraction 
7~W = J2j & ' ) °f components of a vector operator O = (O 1 ,0 2 ,0 3 ) that 

must be traceless, ^ Q J aa = 0, in order to be orthogonal to T^°\ The only available 
vector is the generator of rotations: the spin operator S itself. Its components have 
zero trace because they generate rotation matrices of unit determinant (1 = det U = 
det expli^S*- 7 } = exp{i6* tr 5 J }). Alternatively, one can use the Wigner-Eckart theorem 
[T21 (84)] to show explicitly that S q — ^fc^Tq 1 " 1 up to a normalisation constant such 
that 

= Lj2S j (S j \ ■ ) . (43) 

3 

The normalisation constant c s is fixed by requiring T^-'T^ 1 ' — T^- 1 ': Since S 2 = 
s(s + l)l s is the Casimir operator of the irreducible representation of dimension 
d s = 2s + 1, we have 3 tr {S" i ,S J } = s(s + l)(2s + 1)<%. The normalization factor 
therefore is c s = s(s + l)(2s + l)/3. For electrons, S 1 -? — /2 and Ci/ 2 = \ such that 
^~ — 2 <7 ' 7 ( cr "') ' )■ ^s components in the decoupled basis are indeed those of ©. 
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This completes the derivation of projectors for the scalar electronic diffuson 
vertex. For larger spin, higher orders of K have to be considered which essentially 
involves a Gram-Schmidt procedure. In |Appendix A.2| this is done for the photon 
case s = 1. 



3.4- The crossed vertex as an ordinary operator 

We now turn to the diagonalization of the cooperon vertex X that maps an incident 
tensorial ket product \sa) (g> \sf3) onto the final tensorial ket product \s~f) ® \sS). In 
operator form, 

X =Y \ s ls6)X al3nS (sasi3\ . (44) 

Therefore, the crossed vertex X : H s ® H s ^ H s ® H s can be seen as an ordinary 
linear operator or d 2 s x d 2 s matrix. In this respect, we define the trace for X as: 

Tr c X = ^X Q/3 , a/3 . (45) 

a/3 

This trace definition is invariant under the exchange of spin indices 1(28[) . such that 

Tr L £ = Tr c X (46) 

In other words, the partial Liouville conjugation ^1 \sf3)(s5\ i— > \s5)(s[3\ that maps 
the diffuson onto the cooperon vertex preserves their trace, which in turn will permit 
us to derive useful sum rules between eigenvalues in section T4. 31 

We wish to bring the vertex into a least redundant form for iteration by performing 
the usual Clebsch-Gordan decomposition of the argument and image spaces, i.e., a 
suitable basis change that transforms the direct product D^ s ' <g)£)( s ) of two irreducible 
representations acting on H s (g) H s into the direct sum © ... £)( 2s ) of 
irreducible representations D^ K \ K = 0, 1, ... , 2s. The appropriate recoupling route 
now can be mapped out as 



£>W <g> £>(*') (47) 

\ / 

D (L) 

We first change to the spherical basis \Kq) = ^2 a p(ssaj3\Kq)\sa s/3) of the irreducible 
subspace H^ K \ The corresponding vertex components are 

X Kg , KW = (ssaP\Kq)X a0tlS ( SS1 S\K'q') (48) 

a/3, 7(5 

such that X = YjKq,K'q' \^q)XK q , K'q'(K'q'\. Here, X is decomposed over the 
decoupled operator basis ]Kq)(K'q'\. We can therefore define recoupled irreducible 
operators, as in 128(1 . but in the present context with rank L and their (2L + 1) 
components 

Ti L \K,K') = J2(-) K '~ q '( KK '«-<l'\ Lm ) \ K l)( K 'l'\ ■ ( 49 ) 

qq 1 
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The objects thus used can be precisely located on the map (|T7|) . 

\sa)\sp) (sS\( S1 \ 

\Kq) (K'q>\ (50) 

\ S 

tL l \k,k>) 

The Tm\K,K')'s provide the most natural basis set for exploiting rotational 
symmetries in H s (£> H s . Now the crossed vertex X can be decomposed over this 
basis set: 

X=J2 J2 X Lm(K,K')T^(K,K') (51) 

KK' Lm 

where its irreducible components are 

X Lm (K, K') = J2(-) K '~ q ' ( KK 'l -Q'\Lm) X w , q , . (52) 

qq' 

3.5. Scalar crossed vertex: projectors 

The decomposition into ireducible components is especially profitable for a scalar 
vertex since it has only a single non- vanishing component Xqq (K, K) . The CG- 
coefhcients in (14911 and (|52|l then restrict the sum to K = K' . Therefore, the vertex 
connects irreducible sub-spaces H^ K ^ of equal rank and can indeed be written as 

2s 

x = E ^ T(K) ( 53 ) 

K=0 



with eigenvalues \K = Xqq{K)/^2K + 1 and associated tensors 

T W = V2KTT T (0) (K) =J2\ K( l) ( K <ll ( 54 ) 

These are indeed orthogonal projectors T^ K 'T^ K ) = SkK'T^ k ' onto the irreducible 
subspaces H^ K >. Their matrix elements in the decoupled basis of H s ® H s are 

T SU = x^/to^to, ( 55 ) 

Using ^1 (14a)], one shows that 

TrcTW = J2 T $<*0 = 2K+1 ( 56 ) 

a/3 

and that the projectors sum up to the identity, J2k ^ap\s = &ia&06 ■ 

Even the simplest projector on the singlet space H^ ' has seemingly 

complicated matrix elements in the decoupled basis, 

1 

For electrons, one may check by hand that this is indeed equivalent to the much 
nicer formula (|11[) . For spin 1 particles, the contractions (— ) p 8- p _ q of spherical basis 
components become Sij in the cartesian basis, and one gets 

Tji% = \SilSjk (58) 



T { °Ls = ^{-? 8+P -^-^6- Sa . (57) 
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as used in '8, . These heuristic writings in the decoupled basis are much less systematic 
than the exceedingly simple form l|54|) in the spherical basis. Moreover, the matrix 
elements of directly derived from (|55|l would be rather far from the simple form 
<|12fl we wish to justify. To that purpose, section ETfil discusses a general way of deriving 
the crossed projectors T^ K ^ from the ladder (super-)projectors T' K ' and vice versa. 



3. 6. Recoupling of projectors 

The simple exchange rule C a p n s = X a s,-yi3 provides us with a convenient way of linking 
the diffuson and cooperon projectors. The two diagonalization procedures differ by 
the coupling scheme for the two pairs of spin indices and are related by a simple 
recoupling relation. Indeed, the matrix elements l|40|) of the diffuson projectors, 

7$&6 = Y(-y- 5 {ss 7 -6\Kq}(-y-e{ssa-(3\Kq) (59) 



can be derived from the corresponding cooperon projectors (|55[1 with exchanged spin 
indices /3 <-> 5, T^^g = ^2 q {ssa5\Kq){ss^f3\Kq) , with the help of the appropriate 



recoupling relation |I2. (34)]: 

T St* = E K> ) T Xl> T SU = E R >( K > K ') T S • ( 6 °) 

K' K' 

These relations express a mapping between sets of projector matrix elements 
{TgjJfLg} {T^^p} defined by a transformation matrix R s with elements 



R S (K,K') = (2K + 1)S [ S s S s * | 



(61) 



General 6j-symbol symmetry properties and the orthogonality (|18J) imply that the 
matrix R s is real and circular: R s = R s = R" 1 such that deti? s = ±1 with 
deti?i/2 = det R\ = — 1. This transformation conserves the orthogonality of 
projectors: TY h {T^T^} = (2K+1)5 KK , = Tr^T^T^}. 

Putting the transformations (|60|l to work, the scalar projector l|57|) on the singlet 
state K = is predicted to be given by 

T (0) _ H 2S V^ f sK' T (K') ffi2 s 
a * K , 

Using the electron diffuson projectors (JSJ et iJjJJl, one obtains as expected the singlet 
cooperon projector in the form 1111) . 

<o) _ 1 /V(0) t-(i) ^ 1 



t^, iS = — 2 [n^p - = 2 v^tf - 5 « s5 ^ ■ (63) 

Similarly, the projector H12[) onto the triplet space is 



T^Ls = 3 o^Jrfl + * r «lB ) = * (<W<^ + S aS S 07 ) . (64) 



It-(o) , V(i) \ _ 1 

-a/3, 7(5 " J I 2 J ct<5,7/3 ' g 1 aS^f3 ) ~ 2 (Ua " fL 

This completes the derivation of all projectors for the case of spin i. The photon case 
is discussed in | Appendix A| Now the stage is set for the calculation of eigenvalues. 
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4. Calculation of eigenvalues 

4-1. Diffuson eigenvalues Xk 

The scalar spin superopcrator eigenvalues Xk defined through the general 
decomposition l (57j l are at least (2K + l)-fold degenerate. § They can be calculated 
either by projecting the vertex on an arbitrary component q of the respective subspace, 
Xk = (t^ 



CTq K ^, or directly from the vertex matrix elements £ a f3,j8 as 



_ C m (K, K) _ 1 (K) (R -s 

+ 1 2KTT ^ L i> K « (65) 

= WTT ^ {-) s - {ssa-p\Kq)C a0n5 (-)*-* (ss6 -j\K-q). (66) 

Useful information about the possible form of eigenvalues can be gained from this 
direct calculation. To that purpose, consider an arbitrary microscopic spin interaction 
with matrix elements V ia = (s"f\V\sa) . This interaction can itself be developed in the 
basis of irreducible operators V = J^Kq ^K q Tq K "' '. Its components 

V Kq = (T^ K) \v) = tr {T^V} = J2(-y- a (ssi-a\Kq) V- (a (67) 

are the coupling amplitudes of scalar, vector, quadrupolar type, etc. These amplitudes 
may depend on microscopic degrees of freedom of the scattering object such as the 
orientation of a magnetic impurity. The scattering vertex C a p n s = (V ia (V^)^s) is 
the partial trace over these degrees of freedom, 

C aPn s= (VK q VK^) av n 2s - a - p (s S1 -a\Kq)( S sf3-S\K'-q') (68) 

KqK'q' 

As a scalar vertex has no angular dependence, it always take the generic form 

(VKqV K 'q') av = 0~KK'5qq< S K (69) 

where sk is the vertex eigenvalue in the invariant subspace of rank K for the "vertical" 
coupling scheme. The eigenvalue definition l|66|) then becomes a sum over products of 
four CG-coefficients that defines a 6j-symbol P2J (32)] and finally yields: 

X K = (-) 2s+K EH""' Rs{K', K) s k > ■ (70) 

K' 

Remarkably enough, the only information about the microscopic interaction is carried 
by the coupling constants sk- The 6j-symbol merely provides the recoupling from the 
"vertical" form (cry) «-» (/3S) of the initial product of amplitudes \V ia {V^)ffs) to the 
"horizontal" form (ct(3) <-> (7^) necessary for the diagonalization with respect to the 
multiple scattering product rule (0. 

The simplest possible example of a scalar diffuson vertex is of course given by 
isotropic scattering. The interaction is characterized by Vya — vod^a with vq being 
a complex number. The only effective amplitude is then sq = (2s + l)|«o| 2 , and the 
eigenvalues are simply Xk — \vq\ 2 for all K which is evident at once since this vertex 
is simply proportional to the identity. 

§ Of course, larger degeneracies occur if the vertex possesses even higher symmetries; an elementary 
example is the identity / Q( g i7 5 = S-yaSsp with its single d^-fold degenerate eigenvalue A = 1. 
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The first non-trivial example is given by a vector coupling of the spin-flip form 
(J2J). Using the chosen normalization J3J and S q — \J s(s + l)<i s /3 Tf~' by virtue of 
the Wigner-Eckart theorem |12l (84)], one finds that the only nonzero coupling is the 
vector contribution s% — d s /3. Using l|17f) . the eigenvalues (|66|l are then 

(-)<-+" d < ' ' K\_ K(K + 1) 
( ' s 1 s s 1 J 2s(s + l) ' [U) 

The scalar eigenvalue is identically Ao = 1, for any value of s, as required by trace 
preservation. For electrons (s = ^), we recover moreover the eigenvalue Ai — — | of 
the vector mode K = 1 found heuristically in section I2"3l 

4-2. Cooperon eigenvalues \k 

From the vertex diagonalization 153|) . we find the eigenvalues 

x oo( K ) 1 \ - . . 

XK V^TT zkTT ^ (72) 

= 2j ^ 1 + 1 X! ( ssa / 5 l^>^a/3,7'5( SS 7 , 5|^'7) (73) 

As a function of the elementary coupling coefficients sk defined in l|69|) , the eigenvalues 
are 

XK = (-f s+K £ R *( K '> K ) s «' ■ ( 74 ) 

K' 

Using (|17|) , the eigenvalues of the normalized spin-flip crossed vertex (0} with the only 
non-zero coupling s\ = d s /3 are thus 

*.-(-)-<* +1 >{: : f ( 75 , 

In the electronic case (s = this yields indeed the previously found values Xo = ~1 
in the singlet channel and %i = 5 m the triplet channel. 

^.5. Direct recoupling of eigenvalues 

The precise relation between diffuson and cooperon eigenvalues can be understood by 
observing that both of them are obtained by a recoupling procedure from a "vertical" 
coupling scheme of the initial scattering amplitudes towards the relevant direction of 
diagonalization, "horizontal" for the diffuson vertex, and "diagonal" for the cooperon 
vertex. This implies that the different eigenvalues are linked by simple recoupling 
relations and useful sum rules. 

Starting from a scalar intensity vertex V with no angular dependence, 
characterized by the product {VnqVK'q') av = SKK'Sqq' sk , the interaction vertices 
describing respectively the ladder and the crossed diagrams can be written in the 
form 

£q/3, 7 .5 = SK r sp)i^ X af}jlS = £ s K T^ a . (76) 

K K 

Clearly, the coupling constants sk appear as the vertex eigenvalues for the vertical 
coupling scheme. However, the projectors are not in a form suitable for iteration 
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of the multiple scattering sequence with the horizontal product (J5J. Using the 
transformations 

T *f,l = (-)"E(-) 2S+ "' MK,K')T^ S (77) 

K' 

T Sl = E(-) 2S+K ' *>{K, K>) T$2 S , (78) 

A' 

the vertices can be brought into the suitable form, 

C a/ 3,j8 = E Xk T af}fS X a i3. lS = ^XK . (79) 

K A 

Using the definitions ((77(1 and (|78|l , the eigenvalues A a and \k are then immediately 
derived as function of the vertical eigenvalues sk, given by expressions lf70|l and l(7^|) . 
But this in turn implies that the eigenvalues Xk and \k are also directly linked to each 
other by a simple recoupling procedure. Inverting the relations 1(70(1 and 1(741) with the 
help of the orthogonality relation (|18|l for R s , one obtains the vertical eigenvalues as 

s K = (-) 2s+K EH* 'Rs{k' ; k)x k (so) 

A" 

= J2Rs(K',K) Xk > . (81) 

A" 

Injecting pS jl in fTD ) and (JHU into dl, we indeed find 

X K = Y,R* (K', K) xk' , XK=Y, R * ( K '> K ) X «' ■ ( 82 ) 

K' A' 

These relations had been derived previously in the case of photon scattering (formula 
(52) of 8 , s = 1) and are here generalized to arbitrary spin. These recoupling relations 
replace the heuristic prescription u>2 <-> u>3 for the exchange of contraction weights in 
the photonic case to arbitrary spin. 

Taking the trace l|26(l or (|45|l of the decompositions l|79|l and (|76Jl . we find the 
following useful sum rule: 

5>X + 1)xk = E( 2 ^ + !) A a- = E( 2K + l ) T f3pL = ( 2s + !) s o • (83) 

AAA 

The last equality is explained by the fact that the "horizontal" trace over all modes 
K in (|79|l projects onto the scalar component K = in (|76() . By symmetry of our 
recoupling relations, naturally also the inverse relation holds: taking the "vertical" 
trace in (|7()|l yields 

Y,( 2K + *W = £« a ,p0 = (2s + 1)A . (84) 

K 

Let us demonstrate the power of these relations by taking again the spin- flip vertex 
as a paradigmatic example. Its diffuson and cooperon eigenvalues (|7H and f?5*)l have 
turned out to be equal but of opposite sign, xk(s) = —Xk(s) for any K and s. Indeed, 
comparing the eigenvalue expressions as a function of the elementary couplings s^-, 
l|70ll and l|74() , we can trace back this sign to the fact that the spin- flip vertex has only 
one finite component, the vector coupling s\ = d s /3. The relation 1)82(1 then reduces 
to an orthogonality of 6j-symbols and yields immediately Xk{s) = —Xk(s). The trace 
i|83l reduces to A + 3Ai = (remember s = 0) which immediately fixes the triplet 
eigenvalue to Ai = — | once the trace-preserving eigenvalue Ao = 1 is known. By now 
it should be evident that these values are after all of purely geometrical origin. 
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4-4 ■ Reciprocity 

Quite generally, the diffuson vertex of some microscopic spin interaction V is 
£ Q /3, 7 «5 = (VyaVsp) — (Vy a (V^)ps) while the corresponding crossed vertex is 
X aj 3^s = Ca&^p = (V ia (V^)sf3) av = (V-radV^ps)^- Clearly, these vertices are 
identical if the microscopic interaction is symmetric, V — V * . This complies with the 
general rule that the reciprocity theorem assures perfect equality of ladder and crossed 
contributions if the system's S'-matrix is symmetric [2]. In terms of the irreducible 
amplitudes Vkq defined in l|t)7|) , symmetry of the microscopic interaction is equivalently 
stated as 

V = V* V Kq = (-)*V K - q (85) 
This is to be contrasted with the hermiticity condition, 

V = V i W q = (-) q V K - q (86) 

For example, a simple scalar interaction V^ -* = vqI is hermitian if vq is real. 
Non-hermitian interaction would describe absorption or gain which is known to 
preserve equality between ladder and crossed contributions ■ Indeed, is always 
symmetric since = <5r-o<5<jO\/2s + 1 vq fulfills (55} . 

Less trivially, the spin-flip vector coupling — gJ-S has irreducible 

components 

V$l=5 Kl -Z={-YJ- 9 . (87) 

V c s 

Of course, the interaction is hermitian for a real g. But it is not symmetric in 
general because l|S7)) is different from (— ) q V^_ q as soon as J q ^ (— ) q J- q , i.e., if the 
impurities are not all aligned. This is at the origin of the fact that spin-flip scattering 
suppresses the diffusive pole of the cooperon (the "unitary case" in magnetoresistance 
|5]). Furthermore, it explains the observation that it is precisely the antisymmetric 
part 4^ of the atomic photon scattering tensor (i.e., the vector coupling component) 
that breaks the equivalence of ladder and crossed vertices (20] ■ If, however, one can 
align all vectorial scatterers, for example with an external magnetic field, one can 
choose the quantization axis in this direction such that J q = 5 q oJo, and equality of 
ladder and crossed terms is reestablished. This effect has been observed with electronic 
transport in metal samples containing magnetic impurities and subject to a strong 
magnetic field 21 . There, weak localization corrections to transport are suppressed 
by spin-flip processes at low fields, but are restored at high fields because all magnetic 
impurities are then aligned with the field. A similar argument is at the origin of the 
observed magnetic field enhancement of coherent backscattering of light by a resonant 
sample of ultracold atoms [22J whose groundstate degeneracies are lifted by a magnetic 
field. 

Having at hand the diffuson and cooperon eigenvalues l|70[l and Ij74(l as functions of 
the elementary coupling coefficients sk, we see immediately that a difference between 
eigenvalues is generated only by coupling amplitudes s^- with odd K = 1,3,... (for 
electrons and photons, only K = 1 is possible). A simple inspection of the microscopic 
interaction vertex of a particular physical impurity type permits to decide whether 
this coupling is of vectorial rank K = 1, therefore breaks the equivalence of ladder and 
crossed structures and eventually leads to an effective dephasing of weak localization 
effects, or whether it is of scalar or symmetric type K = 0, 2 and then does not affect 
localization effects. 
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5. Consequences for relevant transport quantities 

5.1. Diffuson spin transport 

The formalism of irreducible spin representations greatly simplifies the expression of 
spin transport quantities that are of interest in the growing number of "spintronics" 
applications. The probability of quantum diffusion, defined in subsection 12.11 with 
spin degrees of freedom can be decomposed into its irreducible components: 

-. 2 s 

P a(} , lS (q, u) = — (G« 7 G^) av = p K(l, • (88) 

The probability Pk(<Z,^) with spatial and temporal Fourier variables q and u 
can be computed independently in each spin sector K if the elementary scattering 
vertex as well as the average propagation between scatterers has been diagonalised 
appropriately. For the summed ladder series T> = C/(l — GC), each diffuson mode up 
to an overall normalisation reads 

4 d) (g^) = 1 i n + K — \ (89) 

1 — Aa'(1 + iwt — Dq i T) 
where the diffusion approximation for the intensity propagation in the Kubo limit 
UT,q£ <C 1 has been made; the diffusion constant in d dimensions is D = £ 2 /to] in 
terms of the scattering mean-free path I and the mean free time r = Ijv — l/(27rpAo) 
(evaluated from the self-energy in the Born approximation with Ao = £ ttQj17 = 1 
in our normalisation). The probability time dependence for each irreducible sector 
therefore is of the form P/f (g, i) ~ exp[— Dq 2 t — t/rd(K)] with a diffuson spin decay 
time 

r d (K)=r-^- (90) 
1 — Xk 

given directly as function of the vertex eigenvalue Xk- The time-dependent diffuson 
probability of classical scattering behaves as 

P$,t,& *)=E e-^"^ ■ (91) 

K=0 

After injection of a certain spin state \sa) at r, the total classical probability 
for final states with arbitrary spin 7 at arbitrary position r' should be normalised to 
unity, Paa,yy(q = 0,t) = 1. Indeed, the trace over the final spin index 7 in (|91|l 
projects onto the scalar component K = 0, such that the conservation of probability 
requires 1/^(0) = which is indeed the case for a unit scalar eigenvalue Ao = 1. 

The overall probability of retaining the initial spin state, say a = +^ =: +, 

is P ++ . ++ (t) = (r| 4 ) ++ ) + c-*/ ri (T| 1 + ) . ++ ) = |(1 + e"'/ 7 "^ 1 )) which deviates only 
for short times from the equidistribution value i. The degree of spin polarisation 

n(t) = P ++t++ (t) - (t) relaxes on the time scale r d (l) as n(t) — e~ t / r<i ^ 1 ^. 

Consider for the sake of concreteness the case of multiple electronic spin-flip 
scattering with the vertex (J2J. This vertex has a negative eigenvalue Ai = — i 
such that formally the decay time TdiX) becomes negative, which renders the above 
predictions inacceptable. In fact, in a disordered electronic sample, there are scalar 
defects responsible for elastic scattering and momentum relaxation, say with a rate 
7 e = n e VQ depending on the density n e of defects and their interaction strength vq . In 
spin space, these vertices have unit eigenvalues both for ladder and crossed vertices. In 
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addition there are magnetic impurities, each with a scattering vertex J2J), centered at 
random positions r m with density n m and corresponding momentum scattering rate 
7m — n m9 2 \J{J + + !)■ The inverse of the total scattering rate 7 = j e + j m is 
the mean-free time r = I/7. The effective average spin-flip scattering vertex then has 
the normalised ladder eigenvalues 

A c /= 7e+7mAg =l-?^(l- \ K ) (92) 
7 7 

which for the spin-flip case are A^ ff = 1 and Af ff = 1 — 47 TO /37. The characteristic 
spin polarisation decay rate reads 

1 _ 1-Af _1-A t lm 4 
r„(l)- r\f ~ r m + ° [ 7 } ~ 3r m ' ^ 

5.2. Cooperon dephasing 

The cooperon is the sum of all maximally crossed diagrams which, strictly speaking, 
starts with the second-order scattering term because the single scattering event is 
already counted in the diffuson: C = XGX/(1 — GX). Here, X is the crossed 
vertex associated to the diffuson vertex for scattering with a total rate 7 = 1/r and 
normalised to a unit diffuson eigenvalue A = 1. Summing the geometric series with 
the returned advanced propagator line (see figure ^c)) gives the decomposition 

2s 

C a p^ s {Q, w)=Y J Ck{Q, u)T^ s (94) 

where Q — k + k' is the sum of external momenta. The cooperon eigenfunctions for 
each irreducible mode are 

t —vjj + DQ Z + 1/t c (K) 
with characteristic cooperon dephasing times 

r c (K) = T-*^ (96) 

1 - XK 

that depend directly on the crossed eigenvalues Xk ■ 

In applications to weak localization, one needs the integrated cooperon, both 
over momenta and spin indices, that counts arbitrary loops of counter-propagating 
amplitudes. Up to a normalisation, 

2s 

P c (w)=^^C a7)7a (Q,u,)= Y, w kJ2 C k(Q^) (97) 

Q 7 K=0 Q 

where the cooperon spin-channel weights wk are determined by the crossed contraction 
wk — Y]~ T^j}ya that corresponds to the sum over all final spin states (3 — 7 in the 
maximally crossed diagrams of figure 1(b). Using the recoupling formula (|16|l . one 
finds 

w K =Y,Rs(K,K>)Y,T^l = Rs(K,0) = (-) 2 ^|^±i . (98) 
K' 7 

Remarkably, although the integrated cooperon describes a renormalisation of intensity 
diffusion (the scalar diffuson mode K — 0), in general all the cooperon modes 
K = 0, ... ,2s contribute with non-zero weights |2"5] . 
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Taking again the case of electronic elastic and spin-flip scattering as an example, 
the singlet and triplet channel weights are Wo = — \ and W\ = §, and the effective 
eigenvalues are 



Xf = 7e+7? " Xg =!-—(! -Xk) • (99) 
7 7 

For a small spin- flip scattering rate 7m/7 ^ 1> t ne corresponding dephasing rates are 

■ <100) 

Inserting the spin- flip eigenvalues xo — — 1 an d Xi = 3> the singlet and triplet 
dephasing times are t c (0) = r m /2 and t c (1) = 3r m /2. If the eigenvalues (|99|l in 
the numerator of l|97|l are approximated by unity, the integrated spin-flip cooperon 
finally reads 

"3 1 11 



2 -iu + DQ 2 + 2/3r m 2 -iuj + DQ 2 + 2/r„ 



.(101) 



The coupling to the uncontrolled degrees of freedom of magnetic impurites dephases 
both the singlet and triplet channels irreversibly and leads to a drastic decrease of 
weak localization effects in disordered electronic samples. A similar effect is found 
for the weak localization of photons scattered by cold atoms with degenerate dipolc 
transitions 

Note that in the case of electronic spin-orbit scattering, the triplet channel with 
positive weight w\ = § is also rapidly damped, whereas the singlet channel with 
negative weight wo = (— ) 2s /(2s + 1) = — ^ survives and leads to antilocalization 
and a positive magnetoresistance 1 1 0| which appears here as characteristic for any 
half-integer spin. 



6. Summary and conclusion 

In this paper, we have developed a systematic method to diagonalize the elementary 
spin scattering vertices which are the building block of diffuson and cooperon multiple 
scattering sequences for particles of arbitrary spin s. Our results therefore provide the 
conceptual background for a truly unified description of the mesoscopic spin physics 
of electrons and photons. We have identified the relevant projectors onto invariant 
subspaces that are irreducible with respect to the rotation group. Once these operators 
have been obtained, the diagonalization allows us to transform the vertical coupling 
scheme for the scattering amplitudes into a horizontal scheme necessary for subsequent 
iteration of the multiple scattering sequence. We have obtained the diffuson and 
cooperon scattering eigenvalues as a function of the microscopic scattering mechanism, 
together with simple recoupling relations as well as useful sum rules. We have shown 
how these eigenvalues directly enter the expressions of the phase coherence times of 
weak localization. 

The method presented here may be extended to non-scalar vertices such as 
the transverse photon propagator G(q) that was diagonalized exactly in H|. This 
transverse propagator is no longer purely scalar at non-zero momentum transfer 
q 7^ 0, but contains quadrupolar parts coupling the modes K = 0, 2. A treatment 
of these non-scalar vertices would start from the expressions (|36|) and (|51|l of the 
present work and derive the appropriate projectors and eigenvalues. On a similar 
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line of thought, light scattering by nematic crystals [23 as well as photon scattering 
vertices of atoms under the influence of an external magnetic field |22| require 
anisotropic diagonalization. Finally, these techniques may become useful for quantum 
computation schemes involving spin degrees of freedom (such as the one studied by 
Loss and di Vincenzo [251 > to cite a paper employing a superoperator formalism quite 
similar to ours) or for entanglement characterization in irreducible representations of 
observable-induced tensor product spaces [26| . 
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Appendix A. Photon scattering 

Appendix A.l. Atomic vertex eigenvalues 

For photons, polarisation-dependent scattering proves more complicated than for 
electrons because of field transversality. The spin degrees of freedom are not decoupled 
from average propagation such that the complete diagonalization of the diffuson and 
cooperon series is much more involved. Nonetheless, the projection onto irreducible 
subspaces permits to derive all eigenvalues and projectors for isotropic photon vertices 
as well, as has been done in J5| for the case of photon scattering from degenerate atomic 
dipolc transitions. For resonant photon scattering, the elementary interaction is of the 
form Vdip = D E where the electric field operator E, proportional to the polarisation 
vector e, creates or annihilates photons whereas the dipole operator D induces internal 
transitions between electronic states with angular momentum J g and J e . Obviously, 
the elementary interaction Vdi p is of the same vectorial type as the electronic spin-flip 
vertex |JT|I. The full photon scattering process, however, comprises the annihilation of 
the incident photon followed by the creation of the scattered one. Therefore, the full 
scattering amplitude is a second-order process in Vdip, and its elementary coupling 
coefficients sk (see eq. (30) of 8 ) contain already all orders K — 0,1,2 obtained 
by the recoupling of two vector interactions of rank one. The ladder and crossed 
eigenvalues are then expressed in terms of 6j- and 9j-symbols. All relations between 
eigenvalues derived in section 0] apply to the photon case as well. The sum rule (|83|l 
has not been evaluated so far. For resonant photon scattering from a closed dipole 
transition J g — > J c , it reads 



The more general case of photon scattering from entire multiplets of hyperfine 
or fine structure dipole transitions can be treated as well |27| . In that case, the 
eigenvalues become frequency dependent, but all algebraic relations of the present 
work continue to hold. 




(A.l) 
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Appendix A. 2. Photon projectors 

Concerning the projectors onto irreducible subspaces, it was shown in [5] that one and 
the same set of projectors 

r^ jk = \sa8 jk , (A.2) 
?u% = l(8ijfa-SikSji), ( A -3) 

"4% = \ (M« + Mil) - I SuSjk. (A.4) 

diagonalises both ladder and crossed vertex alike. Here, the indices i,j,k,l are 
Cartesian indices in configuration space R 3 . 

We still need to connect these tensors to the diffuson (super-)projectors T^ K > 
derived in 13. 3. II and the seemingly different cooperon projectors T' x ' derived in 13.51 
Of course, all formulae of the present work are designed to apply to spin 1 as well. But 
spin 1 is special because its states transform under the fundamental representation 
SO (3) itself. This representation has the peculiar feature that its dimension is 3 and 
therefore equal to the dimension of the abstract group. For spin 1, one therefore can 
use either the standard spherical basis {|lg)} where S z \lq) = q\lq) is diagonal (and 
which was used throughout the present paper) or the Cartesian basis of R 3 itself where 
the generators are S 3 kl = —itijk (in this so-called adjoint representation, the generators 
are essentially given by the structure constants fjki — —£jki of the Lie algebra so(3)). 
Using the latter representation in the definitions (|42|l and l|43|) . we indeed recover 
immediately (|A.2jl and (|A.3|) . 

In order to obtain the symmetric traceless projector K = 2, we have to push the 
calculation one step further. Extrapolating from the cases K = 0, 1, the construction 
rule for higher-rank projectors H38|l should be clear: they are complete contractions 

■j-{K) _ Qhi2-i K ( iii 2 ...«| . ^ (A 5) 

iii 2 ...i K 

of operators O tll2 ' lK that are direct products of K copies of S l with the appropriate 
symmetrisation. Up to normalisation, the operator of the K = 2 projector is 

Q ij = S i S j _ (o fc | SI'S') O k - (0°\ S'S 3 ) O (A.6) 

k 

where \fd~ s O Q = 1 pertains to the scalar projector (|42|l . This construction is nothing 
but the Gram-Schmidt procedure used to orthogonalize the basis {e n } of a vector 
space according to ei i— > e\, e^ i— > ei — [e\ ■ e-i)e\, etc. After short algebra, one finds 

Qi j = 1 ( S i S j + S j S i^ _ s(s + l) ^ i (A.7) 

£ o 
without need for further normalisation. This operator is manifestly symmetric and 
traceless by construction. Using Sti = —i^ijky we obtain the traceless symmetric 
projector l|A.4|) we sought to justify. 

For particles of larger spin than s = 1 (or more general vertices with higher-order 
irreducible components), the calculation of higher-order projectors requires to compute 
traces of increasingly large products of spin operators, tr{S l S : > . . . S p }. Beyond the 
first terms K = 0, 1, 2, carrying out the traces becomes rapidly cumbersome, and the 
formulation (|40|) in terms of spherical components Tq turns out to be more economic 
since the CG-coefficients automatically incorporate the correct symmetrisation. 
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At this point, we have completely justified the diffuson tensors (IA.2[l - i|A.4|) . 
However, at first sight the corresponding crossed projectors T^ K > derived in 13.51 have 
nothing in common with them. But inserting the relevant transformation coefficients 
for s = 1, 



1/3 


-1/3 


1/3 


-1 


1/2 


1/2 


5/3 


5/6 


1/6 



{Ri)k,k> = ( -1 1/2 1/2 (A.8) 

into the recoupling relation l|60(l yields indeed T^j, = T^l, K = 0, 1, 2. This proves 
that the diagonalization of spin 1 ladder and crossed vertices in Cartesian components 
involves the unique set of isotropic projectors QA.2MA.4jl . This is no longer the case 



for half- integer spins s — h, §,•■■ because their SU(2) representations are complex 
unitary which means that ladder and crossed vertices couple rotationally different 
objects. The same conclusion holds for integer spins s = 2,4, . . ., though they admit 
real orthogonal representations of SO (3). This is because the adjoint representation 
is not available and one has to work a priori with two distinct sets of projectors. 
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